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19.  ABSTRACT 


The  equations  of  motion,  in  three  dimensions,  of  solid  particles 
entrained  by  a  gas  flow  through  the  stationary  and  rotating  cascades 
of  a  turbine  are  derived.  'The  gas  velocity  components  and  density  at 
all  the  mesh  points  of  a  square  grid  constructed  in  the  flow  channels 
are  computed  assuming  a  compressible  flow.  Formulas  to  determine  the 
proper  drag  on  the  particles  for  a  wide  range  of  Reynolds  numbers  are 
given.  A  gas  particle  flow  tunnel  ds  used  to  investigate  experimentally 
the  phenomenon  of  particle  impact  with  the  turbine  blades  or  casing  and 
their  rebound  from  these  walls.  Formulas  for  the  restitution  ratio  due 
to  collision  and  the  rebound  to  incidence  angle  ratio  are  derived.  This 
information  is  used  in  the  equations  of  motion  of  the  solid  particles  .7/ 

The  dynamic  behavior  of  the  solid  particles  ip  the  turbine  stage,  i 
namely  their  absolute  and  relative  trajectories,  absolute  nondimens ional 
£  velocity  history  in  the  channel,  and  their  velocity  diagram  as  compared 
to  that  of  the  gas,  is  investigated. 
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ABSTRACT 


The  equations  of  motion,  in  three  dimensions,  of  solid  particles 
entrained  by  a  gas  flow  through  the  stationary  and  rotating  cascades 
of  a  turbine  are  derived.  The  gas  velocity  components  and  density  at 
all  the  mesh  points  of  a  square  grid  constructed  in  the  flow  channels 
are  computed  assuming  a  compressible  flow.  Formulas  to  determine  the 
proper  drag  on  the  particles  for  a  wide  range  of  Reynolds  numbers  are 
given.  A  gas  particle  flow  tunnel  is  used  to  investigate  experimentally 
the  phenomenon  of  particle  impact  with  the  turbine  blades  or  casing  and 
their  rebound  from  these  walls.  Formulas  for  the  restitution  ratio  due 
to  collision  and  che  rebound  to  incidence  angle  ratio  are  derived.  This 
information  is  used  in  the  equations  of  motion  of  the  solid  particles. 

The  dynamic  behavior  of  the  solid  particles  in  the  turbine  stage, 
namely  their  absolute  and  relative  trajectories,  absolute  nondimens ional 
velocity  history  in  the  channel,  and  their  velocity  diagrams  as  compared 
to  that  of  the  gas,  is  investigated.  The  effect  of  different  flow 
parameters,  mainly,  the  particle  mean  diameter,  material  density,  and 
particle  and  gas  initial  velocities  on  the  dynamic  characteristics  of 
the  solid  particles  are  studied. 

Observations  concerning  the  erosion  damage  suffered  by  the  turbine 
stator  and  rotor  blades  as  well  as  the  turbine  casing  due  to  the  solid 
particle  inpingements  are  presented. 
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INTRODUCTION 


Rockets,  aircraft  engines,  and  industrial  gas  turbines  operating 
in  desert  or  dusty  areas  and  in  places  where  the  atmosphere  is  polluted 
by  small  solid  particles  from  factories  and  car  exhausts  can  be  examples 
of  machines  operating  under  gas  particle  two  phase  flow  conditions. 
Another  example  is  engines  that  burn  fuels  which  produce  solid  particles 
in  the  products  of  combustion.  The  solid  particles  mixed  with  the  in¬ 
let  air  or  combustion  gases,  due  to  the  difference  in  their  inertia, 
will  be  driven  away  from  the  streamlines  of  the  gas  and  impact  with  the 
surrounding  walls  of  the  engine.  The  presence  of  solid  particles  in  the 
gas  stream  constitute  a  possible  cause  for  severe  erosion  damage  to  the 
engine  parts.  The  erosion  damage  caused  by  solid  particles  is  expected 
to  be  more  severe  in  the  rotating  parts  of  the  engine,  where  the  particle 
velocities  and  frequency  of  impacts,  as  well  as  the  flow  temperature, 
are  higher.  Hence,  turbines,  axial  or  centrifugal,  of  rockets  or  gas 
turbines  are  parts  of  the  engine  critically  subjected  to  solid  particle 
erosion.  For  example,  in  some  instances  the  life  of  a  turbine  of  a 
helicopter  engine  is  decreased  by  one-fourth  due  to  its  operation  in 
a  solid  particle  suspension  as  compared  to  engines  operating  under 
normal  conditions.  This  gives  an  indication  to  the  seriousness  of  the 
problem  from  the  economical  and  reliability  points  of  view.  The  rate 
of  erosion  damage  to  the  blades  of  the  stator  or  the  rotor  of  turbines 
is  a  function  of  the  blade  and  particle  materials,  the  gas  conditions 
in  the  channel,  the  angle  of  impact  of  the  solid  particles  with  the 
blade  surface,  the  velocity  of  impact  and  the  frequency  of  collisions. 

In  rotating  turbomachinery  the  centrifugal  forces  acting  on  the  par¬ 
ticles,  due  to  its  motion  in  a  circular  path,  tends  to  force  the  solid 
particles  to  move  radially,  and  hence,  impact  with  the  turbine  casing, 
causing  it  to  also  suffer  from  erosion  damage.  Furthermore,  the  rebound 
of  the  particles  from  the  casing  wall  will  increase  the  frequency  of 
collisions  with  the  blades.  This  discussion  points  out  the  fact  that 
the  solid  particles  motion  in  a  rotating  turbomachine  is  a  three- 
dimensional  problem. 

In  order  to  further  understand  the  erosion  phenomenon  of  rotating 
turbomachines ,  it  is  important  to  study  the  dynamic  characteristics  or 
behavior  of  the  solid  particles  entrained  by  the  gas  flow  through  the 
stages  of  turbines.  By  the  dynamic  behavior,  it  is  meant,  the  absolute 
and  relative  position  of  the  solid  particles  everywhere  in  the  channels, 
the  velocity  history  and  velocity  diagrams  of  the  particles  as  well  as 
a  description  of  the  collision  and  rebound  mechanisms  of  the  solid 
particles  from  the  blade  and  casing  walls.  The  particles  trajectories 
studied  in  rotating  turbomachinery  furnishes  information  for  correlating 
experimentally  determined  erosion  severity  of  blades  and  the  character¬ 
istics  of  particle  impact  on  surfaces.  It  also  helps  to  predict  the 
erosion  of  proposed  cascade  designs  based  on  the  correlation  results 
and  thus  initiate  means  to  protect  blades  and  walls  by  minimizing 
erosion  damage.  In  fact,  the  investigation  of  the  solid  particles 
dynamic  behavior  will  enable  engineers  to  design  turbines  and  com¬ 
pressors  to  minimize  erosion  as  well  as  to  optimize  aerodynamic 
characteristics . 

This  investigation  was  devoted  to  studying  the  dynamic  behavior  of 
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solid  particles  in  a  turbine  stage.  The  effect  of  different  flow 
parameters  such  as  particle  mean  diameter,  material  density  and 
initial  particle  and  gas  velocities  on  the  particle  dynamic  behavior 
are  investigated.  Results  of  this  study  are  then  used  to, make  ob¬ 
servations  concerning  the  areas  of  the  blades  or  casing  that  are 
subjected  to  more  impacts  and  hence,  higher  rates  of  erosion  damage. 

I 

The  three-dimensional  equations  of  motion  of  solid  particles 
moving  in  a  compressible  gas  stream  in  a  rotating  cascade  of  a 
turbomachine  are  derived  in  a  general  form.  They  are  solved  to  determine 
the  particle  trajectories  an-,  velocities  in  an  akial  or  radial  rotating 
turbomachinery  or  any  similar  three  or  two-dimensional  particulate  flow 
problem  as  special  cases.  In  order  to  derive  and  solve  the  equations  of 
motion  in  the  three-dimensional  space  of  solid  particles,  several  as¬ 
sumptions  concerning  the  analysis  have  to  be  made.  Also  bassic  informa¬ 
tion  about  the  impact  and  rebound  phenomenon ,  drag  on  particles ,  and 
gas  properties  has  to  be  investigated. 

It  is  assumed  that  the  particles  enter  the  stator  of  the  turbine 
or  compressor  with  uniform  properties  and  equal  distances  apart.  The 
particles  are  assumed  to  be  spheres  of  constant  average  mean  diameter 
uniform  material  density  and  .small  in  size.  The  forces  acting  on  the 
solid  particles  causing  their  motion  in  the  flpw  are  assumed  to  be 
mainly  the  drag  forces  exerted  on  them  by  the  gas.  It  is  further 
assumed  that  the  presence  of  solid  particles  dees  not  alter  the  geis 
properties  from  that  for  the  case  of  gas  flow  alone  passing  through 
the  same  cascade.  This  assumption  is  more  realistic  for  higher  par¬ 
ticle  material  densities  and  small  particle  concentrations.  The 
particle  concentration  «,  is  defined  by  the  ratio  of  mass  flow  rate,  of 
particle  to  the  rate  of  mass  flow  of  the  mixture.  Higher  material 
densities  implies  smaller  number  of  particles  for  the  same  concentration. 

i  : 

After  particles  impact  with  the  blade  or  casing,  they  suffer  a 
drop  in  their  velocities  and  a  change  in  their  direction.  In  order  to 
study  the  impact  and  rebound  phenomenon,  a  simple  test  facility  is  . 

designed,  where  particles  and  blades  of  the  same  materials  as  in  actual 
operating  conditions  are  used.  Particles  are  injected  into  the  flow 
in  the  test  section  and  photographed  by  a  high-speed  camera.  Analysis 
of  the  high-speed  photography  of  the  collision  phenomenon  gives  the 
restitution  ratio  and  die  ratio  of  the  rebound  to  the  incidence  angle  • 
as  a  function  of  the  incidence  angle.  These  two  ratios,  or  any 
equivalent  ratios,  are  sufficient  to  define  the  mechanism  of  impact 
for  a  certain  particle-wall  material  combination.  From  these  ratios 
the  relative  velocity  and  the  direction  of  motion  of  the  particle 
after  collision  are  computed  and  used  as  the  initial  conditions  to  the 
solution  of  the  governing  equations  of  motion  for  the  particles.' 

The  drag  coefficient  on  spherical  particles  has  to  be  known  for 
the  practical  range  of  the  Reynolds  numbers  based  on  particle  mean 
diameter  and  the  magnitude  of  the  difference  in  gas  and  particle, 
velocities.  Analytical  relations  to  determine  the  drag  coefficient  as 
a  function  of  the  Reynolds  number  may  be  derived  from  the  solution  of 
the  Navier-Stokes  equations ,  which  are  valid  for  Reynolds  numbers  up 
to  4.  Empirical  formulas  that  fit  the  experimental  results  for  the 
drag  coefficient  on  spheres  moving  in  a  stream  of  air  are  derived  up 
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to  the  practical  limit  of  the  Reynolds  number  of  such  types  of  flows. 
These  formulas  permit  accurate  calculations  of  the  drag  forces  acting 
on  the  particle  and  increase  the  range  of  application  of  the  equations 
of  motion  of  solid  particles. 

The  remaining  information  that  has  to  be  known  before  solving 
the  particle  equations  of  motion  are  the  gas  properties,  mainly  its 
density  and  velocity  components  everywhere  in  the  channel.  The  equations 
of  motion  of  a  compressible  gas  moving  in  a  rotating  cascade  are  solved. 

A  square  grid  is  constructed  in  the  channel  and  the  gar  properties  are 
calculated  at  all  mesh  points.  The  gas  properties  anywhere  in  the 
nozzle  could  then  be  interpolated.  The  gas  properties  ar;  calculated 
at  all  mesh  points  in  every  row  of  blades  and  stored  on  computer  magnetic 
tape’ in  groups  of  data  each  representing  a  cascade  row  that  constitutes 
the  stages  of  the  turbine. 


The  equations  of  motion  of  the  particles  are  formulated  with  re¬ 
spect  to  axes  fixed  in  the  blades  at  the  entrance  of  the  cascade  row. 

With' the  impact  phenomenon  well  described,  and  the  drag  coefficient  and 
gas  properties  known,  the  equations  of  motion  of  the  particles  are  solved. 
One  row  of  blades  is  considered  at  a  time,  the  particles  enter  the  row 
with  known  initial  conditions,  hit  the  walls  and  rebound  (several  re¬ 
bounds  may  occur)  until  they  leave  the  nozzle.  The  particles  outlet 
conditions  constitute  the  initial  condition  for  the  successive  row. 

The  equations  of  motion  of  particles  with  respect  to  a  new  frame  of 
axes  using  the  known  initial  conditions  are  solved  using  the  correspond¬ 
ing  gas  conditions  in  the  new  cascade  row  and  so  on.  The  particles 
positions  are  referred  to  the  initial  axes  by  a  simple  transformation. 
Repeated  solution  of  the  equations  of  motion  of  the  particles  for  all 
successive  cascade  rows  gives  the  particles  dynamic  behavior  throughout 
the-  turbine.  The  absolute  trajectories  of  the  particles  show  its  real 
path,  while  the  relative  trajectories  indicate  its  path  relative  to  the 
rotor.  This  information  can  be  used  to  determine  the  erosion  rate, 
knowing  the  momentum  loss  due  to  impacts  and  their  frequency  at  a  certain 


rea. 


Tr.e  particles  vs.ocity  aiagrar-s  as  compare c  zo  zr.e  gas  velocity 


ciagram  of  a  scage  vrouic  illustrate  the  general  behavior  of  the  particles 
and  their  deviations  from  the  gas  at  inlet  and  exit  ports. 


The  study  of  particles  trajectories  in  flow  fields  such  as  nozzles, 
pipes  and  cyclones  are  reported  in  many  articles  (References  1-5) .  How¬ 
ever,  none  of  these  investigations  solve  the  -  ..ree-dimcnsional  problem, 
or  take  the  collision  of  particles  with  the  walls  or  the  rotating  cascade 
type  problem  into  consideration.  In  References  6,  7,  and  8  experimental 
and  theoretical  studies  are  made  by  the  author,  to  study  or  simulate  the 
dynamic  behavior  of  solid  particles.  The  results  of  the  two-dimensional 
experimental  study  of  Reference  6,  the  theoretical  study  of  Reference  8 
and  the  results  of  this  investigation  seem  to  agree  with  the  physical 
observation  of  eroded  cascades  of  tested  and  real  engines.  The  results 
of  this  investigation  provide  an  understanding  of  the  behavior  ^ c  the 
particles  in  the  rotors,  a  contribution  that  was  difficult  to  accomplish 
by  the  simulation  of  Reference  7.  The  computer  program  used  in  the  study 
is  given  by  the  authors  in  Reference  9.  It  can  be  used  with  slight 
modifications  to  calculate  particle  trajectories  and  velocities  in  two 
or  three  dimensions,  stationary  or  rotating  axial  or  radial  turbo¬ 
machinery. 
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FORCES  ON  SPHERICAL  PARTICLES  MOVING  IN  A  GAS  STREAM 


The  forces  that  act  on  solid  spherical  particles  suspended  by  the 
gas  flow  are  the  drag  force,  the  force  due  to  the  pressure  gradient  in 
the  fluid  surrounding  the  particles,  the  force  to  accelerate  the  apparent 
mass  of  the  particle  relative  to  the  fluid,  the  Magnus  force,  and  the 
Basset  force  that  takes  into  account  the  deviation  in  the  flow  pattern 
from  steady  state,  see  Reference  10.  Newton's  law  of  motion  may  be  used 
to  descr..be  the  particle  motion,  where  the  sum  of  the  forces  acting  on 
the  particle  plus  its  inertia  force  vanishes.  The  force  on  a  small 
particle  suspended  in  a  turbulent  flow  is  given  in  Reference  11.  Further 
discussions  of  the  problem  are  reported  in  References  12,  13  and  14.  It 
is  assumed  that  the  particles  are  spherical  in  shape,  so  that  the  drag 
force  may  be  calculated  using  the  Stokes  drag  coefficient,  however,  cor¬ 
rection  to  the  value  of  the  drag  forces  at  higher  Reynolds  number  is 
introduced.  It  is  further  assumed  that  the  particles  are  small  when 
compared  to  the  smallest  wave  length  of  the  turbulence,  hence  the  effect 
of  particle  motion  due  to  shear  flow  is  neglected.  Since  in  the  case 
of  particulate  flow  problems  in  gas  turbines,  the  particle  mean  diameter 
is  small  and  the  density  of  particle  material  is  much  higher  in  magnitude 
than  the  gas  density,  the  force  due  to  pressure  gradient,  the  force  to 
accelerate  the  apparent  mass  of  the  particles  and  the  Basset  force  may 
be  neglected  compared  to  drag  force  (Reference  10) .  The  drag  force  and 
force  due  to  pressure  gradient  acting  on  a  spherical  particle  are  dis¬ 
cussed  herein,  however,  the  equations  of  motion  for  the  particles  are 
solved  taking  only  the  drag  for^e  into  account. 


Drag  Force 

An  analytical  solution  of  the  governing  equations  of  motion  for  a 
gas  flow  around  a  spherical  particle  is  available  only  for  very  slow 
moving  particles,  i.e.  for  small  Reynolds  number. 

To  determine  the  drag  force  acting  on  the  solid  particles,  consider 
the  motion  of  a  spherical  particle  in  a  steady,  slow  moving,  incompressible 
gas  flow.  If  the  inertia  forces  can  be  neglected  compared  to  the  viscous 
forces,  and  for  no  body  forces,  the  Navier-Stokes  equations  become; 

grad  p  =  yg  v2  (Cg)  (1) 

The  continuity  equation  of  the  gas  may  be  written  as 

Div  (Cg)  =  0  (2) 

The  solution  of  Equations  (1)  and  (2)  (References  8  and  15)  yields 
the  total  drag  force  acting  on  the  spherical  particles  for  Reynolds  number 
less  than  0.1  as; 

D  =  3  7i  y  d  C  (3) 

g  p  g 

If  the  particle  is  moving  with  an  absolute  velocity  C„  in  a  stream 
of  gas  moving  with  an  absolute  velocity  Cg,  the  drag  force^on  the  spherical 
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(4) 


particle  may  be  written  as; 


5  =  3  TT  j.  d  (C  -  C  ) 

g  p  g  p 


For  higher  Reynolds  numbers,  Equation  (4)  connot  be  considered 
valid,  a  correction  factor  for  this  equation  which  is  a  function  of 
the  Reynolds  number  was  introduced  (References  8  and  15) . 


The  modified  version  of  Equation  (4)  may  take  the  form, 
D  =  3  it  u  d  (C  -  CJ  g  (Re) 

g  p  g  p  3 


where 


g(Re) 


„  r  Re 
D  24 


(5) 


(6) 


In  Equation  (5)  the  Reynolds  number  dependent  function  g(Re)  has 
different  forms  depending  on  the  range  of  the  Reynolds  number.  It  should 
have  a  value  of  one  for  Reynolds  numbers  less  than  1.0.  Its  value  can 
be  determined  theoretically  from  the  solution  of  the  Navier-Stokes 
equations  for  Reynolds  numbers  up  to  4.  The  experimental  data  for  drag 
coefficients  on  a  sphere  has  to  be  used  to  determine  g(Re)  for  Reynolds 
numbers  greater  than  4  and  within  the  practical  limit  of  Reynolds  number 
for  particulate  gas  flow  in  turbomachines.  Equation  (6)  shows  that  the 
function  g(Re)  can  be  given  as  a  ratio  between  the  drag  coefficient  at 
any  Reynolds  number  and  the  Stokes  drag  coefficient  at  small  Reynolds 
number. 

The  formulas  for  the  drag  coefficient  for  spherical  particles  as 
defined  in  Equation  (6)  are  given  below  for  different  ranges  of  the 
Reynolds  number. 


CD  =  g  0  <  Ra  <  1.0  (7) 

A  solution  of  the  Navier-Stokes  equations  of  motion  of  a  solid 
particle  using  a  higher  order  term  gives , 


C 


D 


24 

Re 


(1  + 


1.0  <  Re  <  4.0 


(8) 


The  experimental  data  given  by  Schlichting  (Reference  16)  for  the 
drag  coefficient  of  spherical  particles  are  fitted  (Reference  17)  with  a 
curve  of  the  type 

y  =  a  x*3  +  a  (9) 

where 


yly2  -  y32 

yi  +  y2  -  2y3 


(10) 
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The  resulting  drag  coefficient  formula  is  given  by 


CD  =  21.9416  Re-0*718  +  0.3240 


4.0  <  Re  <_  2000 

For  higher  Reynolds  number  one  may  write , 

CD  =  0.4  Re  >  2000 


(ID 

(12) 


The  Reynolds  number  in  Equations  (7) ,  (8) ,  (11)  and  (12)  is  based 
on  the  particle  mean  diameter  and  the  magnitude  of  the  particle  velocity 
relative  to  the  gas  or. 


i  d  p  d  p 

Re  =  -E_ 2  |  (V  -  v  )  |  =  -£— a 

ii  I  a  n  1  ii 


I  (Sg  -  5p>  I 


or 


Re  = 


d  p 

-p  g 


x)2  +  (v  -  y)2  +  (w  -  z) 

g  g 


(13) 


(14) 


Figure,  1  gives  the  fitted  drag  curve  as  compared  to  the  experimental 
results  an<^  it  shows  good  agreement. 


Force  Due  to  Pressure  Gradient  in  the  Flow  Around  the  Particle 


In  order  to  evaluate  the  force  acting  on  a  spherical  particle  due 
to  the  pressure  gradient  around  it,  consider  a  spherical  particle  of 
diameter  dp  and  center  0  in  a  flow  field  (Figure  2) .  The  pressure  at  any 
point  on  the  particle  surface  is  p(dp/2,  ^ ,  <J>2 )  and  the  pressure  at  0 

is  pQ.  The  polar  coordinates  r,  <pi  and  <t>2  originate  from  the  sphere's 
center.  The  coordinate  $3.  is  measured  from  the  direction  of  the  gradient 
vector  of  the  pressure  at  the  center  (Vp)0.  The  vector  5  is  a  unit  vector 
from  the  center  of  the  sphere  to  any  point  p(dp/2,  <h,  $2)  .  Also  5^  and 
§2  are  the  unit  vectors  in  the  direction  of  (Vp)Q  and  any  normal  to  it  that 
makes  an  angle  (<j>2  -  4>Q)  from  the  projection  of  e  on  the  plane  =  -n/2. 

The  pressure  p (dp/2,  $2)  at  any  point  on  the  surface  of  the  j 

sphere  is  given  by  *  < 

d  d  I 

P  =  P  (2  '  ♦]/  ^2^  =  PQ  +  *  2^  •; 


d 

=  Po  +  -2  -  [(vP}o  *  i] 


P  =  po  +  2  I (Vp)ol  COS  *1 


The  total  force  acting  on  the  particle,  P,  equals  the  sum  of  these 
pressure  forces  on  the  surface  of  the  sphere,  i.e. 


9  =  f  (p  e)  dA 


where  dA  is  the  increment  of  area  of  the  surface  of  the  sphere  whose 
normal  is  S.  Using  Equation  (15)  ,  Equation  (16)  may  be  rewritten  as 

d3  2ir  ;t  _ 

P  =  ~0P-  ]  (Vp)Q|  /  /  e  sin  cop  <#> ^  (d  (d  2 )  (17) 


The  component  of  P  in  the  e.^  direction  is  P^  given  by 

d3  2ir  7T  , 

P-L  =  -g^  |  ( Vp)  Q |  /  /  sin  <)>1  cos  tfrjL  (d  4^)  (d  <j>2) 

00  ~ 


P1  = 


The  component  of  P  in  any  direction  e,  perpendicular  to  P,  is  i>2 
and  may  be  given  by 


P2  8 


I  <Vp>J  1  cos  (*2  ~  lf,o)  (d  +  9*  * 


•  I  sin  cos  (ji,  (d  <j>.)  =  0 


Since  p2  vanishes  for  any  arbitrary  direction  <22*  the  force  P  is 
then  in  the  direction  of  e^,  i.e.,  the  direction  of  the  vector  (Vp)Q  and 
is  given  by 


P  =  /  (p  e)  dA  = 


m 

r&  <Vp)( 

Pr, 
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where  m  and  (m  /pD)  are  the  mass  and  volume  of  one  particle,  respectively, 
p  P  & 

Using  Euler's  equation  of  motion  for  the  gas,  the  force  P,  acting  on 
the  particle  due  to  the  pressure  gradient  may  be  written  as 

m  p 

P  -  (^)  (vp)  =  m  (3  (C  V}  C  (21) 


EQUATIONS  OF  MOTION  OF  SOLID  PARTICLES  ENTRAINED 
BY  THE  GAS  FLOW  IN  A  ROTATING  CASCADE 


The  turbine  may  be  divided  into  separate  successive  rows,  for 
example,  stator,  rotor,  stator,  etc.  The  particles  are  assumed  to  enter 
the  first  row  with  their  initial  conditions  referred  to  a  system  of  axes 
associated  with  the  first  row.  The  conditions  of  the  particles  at  the 
outlet  of  the  first  row  may  be  transformed  to  another  set  of  axes  as¬ 
sociated  with  the  second  row  of  blades,  stationary  or  rotating,  and 
these  will  be  considered  initial  conditions  for  the  second  cascade  row  and 
so  on,  until  the  particles  leave  the  turbine  with  the  conditions  at  the 
exit  of  the  last  cascade  rcw.  In  order  to  derive  the  equations  of  motion 
of  the  particles  suspended  by  the  gas  flow,  their  motion  through  a  ro¬ 
tating  cascade  row  is  considered,  with  the  case  of  particle  motion  in  a 
stator  blade  row,  following  as  a  special  case.  It  is  convenient  to  de¬ 
fine  systems  of  axes  associated  with  each  cascade  row  to  simplify  the 
computer  programming  problems,  as  well  as  yielding  a  simpler  set  of 
equations  of  motion  for  particles  in  the  rotor  blades. 

For  convenience,  the  following  axes  systems  are  defined.  B  is_a 
frame  fixed  in  an  arbitrary  blade  of  the  rotating  cascade  with  ni,  112 
and  n3  as  a  set  of  nonparallel,  noncoplaner,  right-handed  unit  vectors 
in  the  direction  of  the  coordinate  curves  x,  6  and  z,  respectively,  as 
shown  in  Figure  3  and  4.  The  point  B0  is  the  origin  of  the  frame  B, 
which  is  taken  at  the  intersection  of  the  plane  tangent  to  the  blade  row 
at  the  entrance  and  the  blade  axial  chord  in  the  mid-stream  surface  of 
revolution  of  radius  R  that  passes  through  the  middle  of  the  blade 
height.  The  coordinate  curves  x,  0  and  z  are  in  the  axial  or  meridional 
direction,  the  tangential,  and  the  radial  directions,  respectively. 
Further,  it  is  assumed  that  the  frame  B,  fixed  in  the  blade  row,  moves 
with  a  constant  angular  velocity  co,  equal  to  the  angular  speed  of  the 
rotor,  in  a  reference  frame  E  fixed  in  the  engine.  The  mutually 
perpendicular  set  of  axes  X,  Y  and  Z,  are  fixed  in  the  engine  at  a 
point  E0.  Point  E0  is  the  intersection  of  the  engine  axis  with  the 
plane  tangent  to  the  blade  row  at  the  entrance.  The  coordinate  axis 
X  is  in  the  axial  direction,  while  Y  and  Z  are  axes  normal  to  the  X  axis 
in  the  plane  tangent  to  the  blade  row  at  the  entrance,  as  shown  in 
Figures  3  and  4. 

Consider  the  motion  of  a  solid  particle  entrained  by  the  gas  flow 
past  a  rotating  cascade.  Referring  to  Figure  5,  the  coordinate  curves 
x,  8  and  z,  that  rotate  with  an  angular  velocity  w,  after  a  time  t,  are 
at  an  angle  equal  to  cut  from  the  Z  axis.  The  particle  p  is  at  any 
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arbitrary  position  x,  0  and  z  measured  from  B0. 

According  to  Newton's  law  of  motion  the  forces  acting  on  the 
particle,  taken  to  be  the  drag  for-a,  equals  the  mass  of  the  particle 
multiplied  by  its  acceleration  in  an  inertial  reference  frame.  The 
equations  of  motion  of  the  particle  p  with  respect  to  the  axes  Xf  Y 
and  Z  of  the  fixed  reference  frame  E  are 


Vp  =  DX 

mp^p  =  ^Y 

VP  -  Ds  <22> 

where  mp  is  the  mass  of  one  particle  given  by 

\  -  k  *  d3P  5P  <23> 

and  Xp,  Yp  and  Zp  are  the  second  derivatives  with  respect  to  time  of 
the  components  of  the  position  vector  of  p  in  E. 

From  Equation  (5)  the  components  of  the  drag  force  in  the  reference 
frame  E,  as  given  in  Equation  (22)  are  written  as 

Dx  =  3  *  ug  dp  g (Re)  (Xg  -  Xp) 


Ov  =  3  *  d„  g (Re)  (Y  -  Y  ) 


g  p 


g  p 


Dz  =  3  u  yg  dp  g(Re)  (Zg  -  Zp) 


In  Equation  (24)  Xg,  Yg  and  Zg  are  the  gas  velocity  components  at 
p  in  E,  while  Xp,  STp  and  ip  are  the  velocity  components  of  a  particle  at 
p  in  E, 

Substituting  Equations (53)  and  (24)  into  Equation  (22)  we  get 

*p  =  G<*g  '  V  (25) 


where 


1  i>  y 

G  =  — S-  g  (Re) 

cr  p 
P  P 


The  coefficient  G  is  a  measure  of  the  particle  characteristic  time 
which  indicates  the  relaxation  time  of  the  particles  relative  to  the 
gas  flow. 

From  Figure  5,  the  relation  between  the  components  of  the  position 
vector  and  its  first  and  second  derivatives  with  respect  to  time  measured 
from  Eq  and  the  components  of  the  position  vector  and  its  time  deriva¬ 
tives  measured  from  Bo  of  the  particle  p,  is  given  by 


X  -  x 


X  »  X 


X  =  x 


Y  =  (R  +  z)  sin  (0  +  to t) 

Y  =  z  sin  (0  +  <ot)  +  (R  +  z)  (0  +  <o)  cos  (0  +  tot) 

Y  =  z  sin  (0  +  tot)  +  2z  (0  +  to)  cos  (O  +  tot) 

+  (R  +  z)  t-(6  +  <o)2  sin  (0  +  tot)  +  0  cos  (0  +  tot)) 


=  (R  +  z)  cos  (0  +  <ot) 

=  z  cos  (0  +  tot)  -  (R  +  z)  (0  +  to)  sin(0  +  tot) 

»*  c  • 

=  z  cos  (0  +  tot)  -  2z  (0  +  to)  sin  (0  +  tot) 

-  (R  +  z)  [(0  +  to)2  cos  (0  +  tot)  +  0  sin  (0  +  tot)] 


The  velocity  components  of  the  gas  at  any  point  p  with  reference 
to  E0,  neglecting  the  gas  velocity  in  the  radial  direction  compared 
to  the  gas  velocity  in  the  axial  and  tangential  directions,  are  written 
as, 

K  = 

g  g 

Y  =  [Vg  +  (R  +  z)to]  cos  (0  +  tot) 

Z  =  [v  +  (R  +  z)  to]  sin  (0  +  tot)  (32) 

g  g 


where  Ug  and  vg  are  the  components  of  gas  velocity  relative  to  the 
rotating  blades  in  the  axial  and  tangential  directions,  respectively. 

Substituting  Equations  (29)  ,  (30) ,  (31)  and  (32)  in  Equations  (25) , 
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(26)  and  (27),  find 


x  =  5  (ug  -  x)  (33) 

z  sin  (0  +  wt)  +  2z  (0  +  to)  cos  (0  +  tot) 

-  (R  +  z)  [{0  +  to)^  sin  (0  +  tot) 

+  0  cos  (0  +  tot)  ]  =  G{  [Vg.  +  (R  +  z)io]  cos  (0  +  tot) 

-  z  sin  (0  +  tot)  -  (R  +  z)  (0  +  to)  cos  (0  +  tot)}  (34) 

z  cos  (0  +  wt)  -  2z  (9  +  to)  sin  (0  +  tot)  -  (R  +  z) 

O  •• 

[(0  +  to;  cos  (0  +  tot)  +  0  sin  (0  +  tot)] 

=  G{  [v  +  (R  +  z)to]  sin  (0  +  tot)  -  z  cos  (0  +  tot) 
y 

+  (R  +  z)  (0  +  to)  sin  (6  +  <ot)  }  (35) 


Multiplying  Equation  (34)  by  cos  (0  +  tot).  Equation  (35)  by 
sin  (0  +  tot)  ,  and  .subtracting  the  two  resulting  equations  and  rearrang 
ing,  we  get  the  equation  of  motion  of  the  particle  p  in  the  tangential 
direction  as 


(R  +  z)  0  =  G[v  -  (R  +  z)  0 ]  -  2z  (0  +  to)  .  (36) 

y 

Multiplying  Equation  (34)  by  sin  (8  +  tot).  Equation  (35)  by 
cos  (8  +  tot)  ,  and  adding  the  two  resulting  equations  and  rearranging, 
we  get  the  equation  of  motion  of  the  particle  p  in  the  radial  direction 
as 


z  =  -Gz+(R+z)  (0  +  to)2  .  (37) 

The  distance  y  that  a  particle  travels  in  the  0  direction  on  the 
surface  of  revolution  with  radius  (R  +  z)  is  given  by 

y  =  (R  +  z)  a  (38; 

For  convenience  Equations  (33),  (36)  and  (37)  may  be  rewriteen  as 

x  =  G  (Ug  -  x)  (39) 

8  =  TR-r-zT  tvg  -  (R  +  to  +  »>  (40) 

z  =  -  G  z  (R  +  z)  (§  +  to)2 
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(41) 


Re  and  G  are  given  by  Equations  (14)  and  (28)  ,  respectively. 


Equations  (39)  ,  (40)  and  (41)  are  the  governing  equations  of 
motion  of  any  particle  p  moving  with  the  gas  stream  through  the  nozzles 
of  a  rotating  cascade.  They  represent  the  equations  of  motion  in  the 
axial,  tangential  and  radial  directions,  respectively.  The  coordinate 
curves  x,  6  and  z  measured  in  the  frame  of  axes  B  and  their  first  and 
second  time  derivatives  as  well  as  the  relative  gas  velocity  components 
in  the  axial  and  tangential  directions  are  taken  at  the  particular  point 
under  consideration.  It  is  noted  that  the  second  term  in  the  right  hand 
side  of  Equation  (40)  represents  Coriolis  acceleration,,  while  the  second 
term  in  the  right  hand  side  of  Equation  (41)  represents  the  centrifugal 
accelerations  of  a  particle  moving  in  curved  path  in  the  rotating  frame 
B.  (References  18  and  19) .  These  equations  form  a  system  of  nonlinear 
differential  equations.  They  are  solved  numerically  for  every  particle 
entering  the  cascade  in  increments  oi  time  if  the  initial  conditions  of 
the  particle  are  known.  The  time  increment  may  have  a  constant  value 
as  long  as  it  does  not  take  the  particle  beyond  the  vail,  thus,  at  the 
regular  increment  nearest  the  wall,  the  time  increment  has  to  be  iterated 
to  the  exact  value  that  is  necessary  for  the  particle  to  just  hit  the 
wall.  Solution  of  the  general  system  of  equations  of  the  type  of  Equa¬ 
tions  (39) ,  (40)  and  (41)  and  a  discussion  of  the  error  are  given  in 
Reference  9.  The  solution  of  the  particle  equations  of  motion  gives 
the  relative  and  absolute  locations  as  well  as  the  velocity  components 
of  a  particle  in  the  rotating  cascade  row.  Successive  solution  of  these 
equations  for  every  cascade  row  will  give  the  particles  path  throughout 
the  turbine.  The  solution  of  these  equations  requires  the  knowledge  of 
the  coefficient  G,  the  collision  and  rebound  mechanism  of  the  particles 
and  the  turbine  walls  as  well  as  the  gas  properties  in  the  channel. 


EXPERIMENTAL  STUDY  OF  THE  IMPACT  AND  REBOUND 


PHENOMENON  OF  SOLID  PARTICLES  FROM  THE  WALLS 


When  a  solid  particle  moves  in  a  stream  of  gas,  it  does  not  in  general 
follow  the  streamlines  taken  by  the  gas  due  to  its  higher  inertia.  For 
particulate  gas  flow  passing  through  the  nozzles  of  a  rotating  cascade, 
j  the  particles  tend  to  collide  with  the  blade  surfaces.  They  often  collide 

i  with  the  turbine  casing  due  to  the  effect  of  the  centrifugal  force  act¬ 

ing  on  the  particle.  After  hitting  the  walls,  the  particles  experience 
a  loss  in  their  momentum  relative  to  the  wall  and  change  the  direction 
of  their  motion.  They  may  hit  the  blades  one  or  more  times  before 
leaving  the  cascade.  These  repeated  collisions  of  several  particles 
causes  severe  erosion  damage  to  the  blades.  The  analytical  study  of 
the  impact  phenomenon  of  particles  with  a  rotating  blade  in  three 
dimensions  is  a  very  difficult  task.  The  value  of  the  particle  velocity 
and  direction  of  its  motion  as  it  rebounds  from  the  surface  after 
collision  must  be  known  in  order  that  the  solution  of  the  particle 
equations  of  motion  be  continued  beyond  the  points  of  collision.  In 
order  to  investigate  the  collision  phenomenon,  an  experin  ntal  study 
was  made,  where  the  particles  were  photographed  and  their  behavior 
during  this  process  is  studied.  The  experimental  facilities  of  Reference 
15  were  used  and  extensive  results  for  steel  blades  and  corn  cups  particles 
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have  been  collected  and  analysed  so  as  to  determine  the  behavior  of  solid 
particles  before  and  after  the  collis:  on  with  the  walls.  The  result  of 
this  study  is  given  in  details  in  References  8  and  15  and  may  be  summarized 
as  follows. 

The  change  in  particle  momentum  u.:e  to  collisions  was  found  to 
be  mainly  a  function  of  the  particle  incidence  angle  0i.  In  Figure  6, 
the  ratio  of  the  particle  velocity  tangent  to  the  surface  after  and 
before  collision  (V  /V  )  is  plotted  versus  0,,  and  may  be  considered 

pt2  ptl  x 

constant  for  all  particle  diameters.  The,  dotted  line  in  Figure  6 
represents  the  straight  line  fitting  the  experimental  data  expressed 
by  the  equation. 


V  /V  =  0.95  +  0.00055  0,  (42) 

pt2  '  ptl  1 


The  change  in  the  ratio  betwee  i  particle  velocity  components  normal 
to  the  surface  of  the  blade  after  and  before  collision  (V  /V  ) ,  as  a 

pn2  pnl 

function  of  8i,  is  given  in  Figure  7  for  different  particle  diameters. 
The  dotted  curve  in  Figure  7  represents  the  curve  fitting  of  the  data 
points  on  the  figure  expressed  by  the  following  relation 


1.0  -  0.002108  0 ^  +  0.0001417  ** 


(43) 


The  restitution  ratio  (V  /V  )  is  defined  as  the  ratio  between 

p2  P1 

the  particle  velocities  after  and  before  collision.  Another  important 
factor  is  the  rebound  to  incidence  angle  ratio  (82/81) >  defined  as  the 
ratio  of  the  angle  between  the  direction  of  the  incidence  particle 
velocity  and  the  tangent  to  the  surface  at  the  point  of  impact  and  the 
angle  between  the  direction  of  the  rebound  particle  velocity  and  the 
tangent  to  the  surface  at  the  same  point. 

Figure  8  gives  82/81  as  a  function  of  81.  The  dotted  curve  is 
obtained  by  substituting  Equations  (42'  and  (43)  into  the  following 
relation. 


89  1  pt2 

ir  *rot " 

1  1  pti 


cot  81J 


In  Figure  9,  the  restitution  ratio  is  presented  versus  83..  Equa¬ 
tions  (42)  and  (43)  together  with  the  following  equation  are  used  to 
plot  the  dotted  curve  in  the  figure. 
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Once  the  particles  and  blade  materials  are  known  for  a  certain  gas 
turbine  application,  fundamental  experiments  can  be  conducted  to  render 
information  about  the  impact  and  rebound  phenomenon  of  this  particular 
material  combination.  It  may  be  added  that  in  order  to  simplify  the 
analysis  the  formulas  found  to  describe  the  impact  phenomenon  do  not 
take  into  account  any  deposits  or  wear  that  may  be  suffered  by  the  blade 
material  due  to  collisions.  Either  Equations  (42)  and  (43)  or  Equations 
(44)  and  (45)  are  enough  to  define  the  particle  condition  just  after 
rebound  from  the  surface.  They  give  new  initial  conditions  to  be  used  in 
the  solution  of  the  equations  of  motion  of  the  particles  after  collision. 


GAS  FLOW  PROPERTIES  IN  A  BLADE  TO  BLADE  SURFACE 
OF  REVOLUTION  OF  A  ROTATING  CASCADE 


The  assumption  that  the  gas  flow  moves  in  the  blade  to  blade 
surface  of  revolution  shown  in  Figure  10 ,  is  applicable  to  axial  or 
radial  flow  cascades.  The  coordinate  curves  x,  8  and  z  are  fixed  in 
the  blade  as  described  earlier,  in  general  they  are  in  the  meridional 
direction,  tangential  and  normal  to  the  blade  to  blade  surface  of 
revolution.  It  is  required  to  determine  the  properties  of  the  com¬ 
pressible  gas  flow,  namely  the  gas  density,  pg,  and  the  velocity 
components  in  the  meridional  and  tangential  direction,  Ug  and  vg,  in 
both  the  rotating  channels  of  the  cascade  and  on  the  blades.  A 
solution  to  the  momentum,  continuity,  energy  and  state  equations  that 
conveniently  describe  the  gas  flow  ir.  a  rotating  cascade  provide  the 
gas  properties  required  in  the  solution  of  the  governing  equations  of 
motion  of  the  particles.  For  small  particle  concentrations  it  can  be 
assumed  that  the  gas  conditions  calculated  from  the  solution  of  the 
nonparticulate  gas  flow  equations  of  motion  in  a  rotating  turbomachine 
are  not  going  to  be  altered  due  to  the  presence  of  particles  in  the 
flow. 


For  a  compressible,  nonviscous  and  steady  gas  flow,  the  momentum 
and  continuity  equations  (References  15  and  20)  may  be  written  as. 
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y  is  the  modified  version  of  the  stream  function. 
The  combined  energy  and  state  equation  is 


|C  |2 

31 
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'm 


(48) 


The  turbine  cascade  channel  region  ABCDEFGH  in  Figure  11  is 
divided  into  square  grids,  the  difference  equation  (Reference  9)  of 
Equation  (46)  is  written  at  all  the  mesh  points.  A  solution  by  the 
method  of  over- relaxation  (Reference  21)  to  the  n  resulting  algebraic 
equations  is  obtained  using  assumed  initial  values  for  the  gas  density, 
say  the  absolute  total  density  at  inlet.  The  solution  of  these  equa¬ 
tions  gives  the  modified  stream  function  at  the  mesh  points.  The 
approximate  relative  gas  velocity  components  at  all  the  mesh  points 
are  then  calculated  from  Equations  (47).  Equation  (48)  may  then  be 
used  to  compute  new  gas  densities  at  the  mesh  points.  These  new 
densities  are  used  to  determine  a  better  approximation  for  the  modi¬ 
fied  stream  function  and  hence  the  relative  gas  velocity,  which  are 
then  used  to  determine  new  gas  densities  and  so  on  ,  until  the  required 
accuracy  in  calculating  the  gas  density  is  reached.  The  fortran 
program  Reference  22  is  modified  slightly  to  give  the  relative  gas 
velocity  components  and  densities  at  the  mesh  points,  by  solving 
L^uations  (46)  and  (48)  then  store  the  data  for  a  cascade  row  on 
magnetic  tape,  to  be  used  as  data  in  the  solution  of  the  particle 
equations  of  motion.  Other  methods  to  determine  the  gas  properties 
in  the  cascade  by  placing  singularities  on  the  contour  of  the  airfoil 
(Reference  23) ,  or  by  distributing  the  singularities  on  the  camber 
and  chord  lines  (Reference  24)  may  be  used.  The  gas  flow  in  these 
cases  is  considered  incompressible  and  correction  for  the  compres¬ 
sibility  effect  has  to  be  introduced  (Reference  25) .  Another  method 
(Reference  26)  gives  the  compressible  gas  flow  properties  for  a  sub¬ 
sonic  cascade  flow  by  the  method  of  distributing  vorticities  along  the 
blade  contour.  The  advantage  of  using  a  square  grid  is  that  it  may 
be  used  as  a  means  to  specify  the  particle  location  in  the  channel. 

The  decision  on  which  method  to  use  to  determine  the  gas  properties 
should  be  a  compromise  between  accuracy,  range  of  application  and 
computer  time. 


NUMERICAL  EXAMPLE 

The  dynamic  characteristics  of  particles  entrained  by  a  gas 
through  a  turbine  stage  consisting  of  a  stator  followed  by  a  rotor 
is  considered  herein.  The  effect  of  different  particle  and  gas  flow 
parameters  on  the  particle  characteristics  are  investigated.  The 
parameters  considered  are  the  particles  mean  diameter,  their  material 
density  and  their  initial  nondimens ional  absolute  velocity,  as  well  as 
the  gas  velocity  at  the  inlet.  The  particle  trajectories  and  velocities 
are  calculated  from  Equations  (38) ,  (40)  and  (41) .  The  drag  coefficient 
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is  calculated  using  Equations  (7),  (8),  (11)  and  (12),  the  impact  and 
rebound  characteristics  from  Equations  (44)  and  (45) ,  and  the  gas 
properties  in  the  cascade  are  obtained  by  solving  Equations  (46)  and 
(48) .  The  computer  program  described  in  Reference  9  is  used  on  the 
IBM  360  computer,  and  the  output  punched  on  cards.  These  cards  are 
fed  as  an  input  to  the  Calcomp  plotter  which  is  programmed  through  the 
IBM  1130  computer  to  plot  the  trajectories ■ and  velocities'  of  the 
particles. 


Particles  Dynamic  Behavior  in  the  Turbine  Stage 


The  dimensions  of  che  stator  and  rotor  of  the  turbine  stage  are 
given  in  Figure  12,  while  the  airfoil  dimensions  are  given  in  Reference  27. 
Figure  13  illustrates  the  combined  gas  velocity  diagram  for  the  53.42% 
degree  of  reaction  turbine  stage.  The  mesh  structure  used  for  the 
calculation  of  the  gas  properties  and  the  particle  trajectories  is  shown 
in  Figure  14.  The  particles  'are  assumed  to  enter  the'  turbine  stator 
with  uniform  velocity.  The  particle  concentration  a  does  not  appear 
in  the  equations  of  motion  of  the  particles,  hence  it  does  not  ir*'- 
fluence  the  behavior  of  any  individual  particle  for  practically  small 
concentrations.  Associated  with  certain  a,  and  suspension  flow  rate 
are  specific  npacings  between  particles,  which  determine  the  place  and 
the  frequency  of  collisions  and  consequently  affects  the  rhte  of  blade 
erosion  at  the  different  points  on  its  surface.  Since  the  blade  .  1 

rotational  speed  is  much  higher  than  the  particle  velocity:  at  the 
entrance  to  the  turbine  rotor,  the  position  of  the  rotor  blades  may 
be  taken  arbitrarily  (Figure  12) .  The  dynamic  behavior  of  a  group  : 
of  five  particles  of  mean  diameter  dp  =  200  microns,  entering  the 
stator  blade  uniformly  is  shown  in  Figures  15,  16,  17  ancl  18.  In  ' 

particular,  the  projection  of  the  absolute  and  relative  trajectories 
in  the  stream  surface  of  revolution  and  the  radial  displacement  and  i 

the  nondimensional  absolute  particle  velocity  distribution  through 
the  stage  is  illustrated  in  these  figures.  Similar  figures  for  eleven 
particles  uniformly  entering  the  stator  cascade  with  mean  diameters  40 
microns  and  3  microns  are  given  dn  Figures  19  to  22  and  23  to  PC,  re¬ 


spectively.  l'n  the  case- of  the  particles  of  dp  =  200  microns,  the 


results  are  only  plotted  for  every  other  particle  in  order  to  facilitate 
interpretation  of_the  figures.  In  Figures  15  through  26,  the  particle 
material  density  pp  is  68.7  Ib/ft3  and  the  initial*  particle  absolute 
velocity  ratio,  nondimensionalized  with  respeqt  to  gas  inlet  velocity, 

C  /C  equals  0.3.  The  gas  conditions  at  the  turbine  stator  inlet 
P  -i  9-i  i 


x 
are 


gas  initial  velocity  C  =  142.64  ft/sec,  initial  gas  density 


P  =  0.076  lb/ft  ,  initial  gas  total  temperature  T  =  6Q°F,  and  the 
9-5  9 


in 

total  mass  flow  rate  per  channel  W  =  0.123  Ib/sec. ;  and  the  particles  ' 
enter  at  the  stream  surface  of  revolution  of  radius  7.05  inches  The 
diameter  of  the  turbine  is  taken  relatively  small  so  that  the  centrifugal 
effects  on  the  particles  are  pronounced  and  in  order  to  give  a  more 
general  particle  path.  The  cascade  design  angular  speed  at  the  mid¬ 
radius  R  is  taken  to  be  ua  =  603.5  radians  per  second, 1  which  is  equivalent 
to  a  mid-radius  blade  velocity  of  352.5  ft/sec.  From  Figures'  15  through 
26  it  can  be  observed  that  as  the  particle  mean  diameter  decreases,  the 


particles  experience  less  deviation  from  the  gas  streamlines  in  the 
stator.  They  tend  to  enter  the  opposite  rotor  blade  channel  and 
experience  less  scattering  in  the  rotating  nozzles  (Figures  16,  20 
and  24) .  Particles  with  higher  mean  diameters  hit  the  blunt  leading 
edge  of  the  rotor  blade  and  return  to  the  stator,  collide  with  the 
stator  suction  surface  and  return  to  the  rotor.  They  may  repeat  the 
same  process  more  than  one  time,  causing  erosion  of  the  suction  side 
of  the  stator  blades  (particles  number  2,  4,  6  and  8  in  Figure  15). 

The  particles  enter  the  rotor  when  they  pick  up  enough  speed  to  take 
them  through  the  rotor  nozzles.  The  band  that  contains  the  particles 
as  they  leave  the  stator  narrows  as  the  particle  mean  diameter  de¬ 
creases  (compare  Figures  15  and  23  for  dp  =  200  and  8  microns,  respec¬ 
tively).  As  the  particles  mean  diameter  decrease,  they  attain  higher 
velocities  in  the  stator  and  the  change  in  their  absolute  velocities 
as  they  enter  the  rotor  decrease  (Figures  18,  22  and  26).  Particles 
with  smaller  mean  diameters  have  smoother  absolute  and  relative  paths 
ar.d  velocity  distributions,  as  can  be  seen  from  Figures  15  through  26. 

The  velocity  distributions  of  particles  with  smaller  dp,  however,  are 
'more  affected  by  the  gas  conditions,  particularly  near  the  stagnation 
points  (Figure  26  for  particles  with  dp  =  8  microns) .  Particles  tend 
to 'move  with  a  higher  radial  velocity  after  hitting  the  rotor  blades. 

The  higher  the  va3  ue  of  dp,  the  higher  will  be  the  radial  velocities 
of  the  particles  after  collision  (Figures  17  and  25) .  The  particles 
travel  toward  the  blade  tip  as  they  pass  through  the  stage.  Those 
particles  of  dp  =  200  microns  hit  the  turbine  casing  before  leaving  the 

turbine  stage  (Figure  17).  Particles  with  dp  =  8  microns,  however, 

travel  through  the  stage  without  hitting  the1  casing  (Figure  25).  The 
vertical  lines  in  the  relative  trajectories  of  Figure  16  indicate  the 
distances  traveled  by  the  rotor  blades  during  the  time  taken  by  the 

particles  to  re-enter  the  stator,  hit  the  stator  blade  suction  side 

and  reach  the  stator  exit.  These  distances  differ  with  the  particle 
initial  position  and  velocity  as  it  returns  to  the  stator.  The  loops 
in  the  curves  of  the  particles  velocity  ir.  Figure  18  are  associated 
with  the  paxticJes  that  return  to  the  stator.  These  loops  result  from 
the  particles  momentum  gain  due  to  their  collision  with  the  rotor  blades, 
followed  by  a  momentum  loss  resulting  from  their  impacts  with  the  stator. 
Particles  with  dD  greater  than  200  microns,  tend  to  move  back  and  forth 
many  times  from  totor  to  the  stator  and  vice  versa,  until  they  either 
break  into  smaller  particles,  or  reach  the  proper  speed  and  direction 
to  go  through  the  rotor.  It  is  also  observed  from  Figures  15  and  16 
for  particles  with  dp  =  200  microns  and  Figures  23  and  24  for  particles 
with  dp  =  3  microns  that,  particles  with  smaller  mean  diameters  experience 
fewer  collisions  with  the  stator  and  the  rotor,  particularly  on  the 
suction  sides  of  the  blades,  ’’article  collisions  with  the  blades  due 
to  their  lateral  motion  are  affected  by  the  cascade  geometry,  namely, 
the  pitch  and  blade  camber.  Other  factors  that  influence  these  collisions 
are  the  mechanisms  of  impact  and  rebound,  the  cascade  rotational  speed, 
and  gas  and  particle  properties.  The  parts  of  the  stator  blade  sub¬ 
jected  to  more  impacts,  and  hence,  higher  rates  of  erosion,  are  the 
leading  ^dge  and  pressure  side.  The  stator  blade  suction  side  is 
subjected  to  impacts  only  from  particles  of  higher  mean  diameters  that 
return  to  the  stator  after  hitting  the  rotor  blade  leading  edge  (Figure 
•15).  In  the  rotor  blades,  the  particles  hit  the  blades  with  higher 
relacivo  velocities,  nence ,  the  rate  of  erosion  is  expected  to  be  higher 
theLn  that  for  the  stator  blades.  The  rotor  blade  leading  edge  is  eroded 
by  particles  of  all  sizes,  while  the  pressure  side  is  subjected  to 
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collisions  only  from  particles  with  higher  mean  diameters  (Figures  16, 

20  and  24) .  As  the  particles  reach  the  casing,  they  impact  and  rebound 
in  an  occilating  manner  v/ith  decreasing  amplitude.  Finally,  they 
move  into  the  stream  surface  of  revolution  that  passes  through  the 
blades  tips  and  cause  severe  damage  to  blade  tips  in  the  later  stages. 

The  previous  conclusions  concerning  the  blade  areas  subjected  to  erosion 
are  in  agreement  with  the  observations  of  actual  eroded  turbine  blades. 

In  order  to  study  the  effects  of  different  particulate  flow  para¬ 
meters  on  the  dynamic  behavior  of  the  particles,  the  particles  are 
divided  into  three  typical  trajectory  groups  according  to  their  location 
at  the  time  they  enter  the  stator.  Figures  15,  19  and  23  or  Figures  27, 
32  and  37  illustrate  typical  examples  of  the  three  trajectory  groups. 

In  particular,  particles  such  as  number  2  start  by  hitting  the  upper 
part  of  the  stator  blade  leading  edge,  particle  number  6  moves  straight 
in  the  stator  channel  and  hits  its  pressure  side,  while  particle  number 
10  starts  by  hitting  the  lower  part  of  the  turbine  stator  leading  edge. 

The  dynamic  behavior  of  particles  having  mean  diameters  of  200,  40 
and  8  microns,  the  same  material  density  and  initial  velocities,  is 
investigated.  The  results  are  shown  in  Figures  27  to  41.  The  effect 
of  particle  mean  diameter  on  it  dynamic  behavior  is  discussed  above. 
Additional  figures, showing  the  combined  velocity  diagrams  of  the  par¬ 
ticles  as  compared  to  those  of  the  gas,  are  given  for  all  of  the 
typical  particles  (Figures  21,  36  and  31).  These  figures  illustrate  the 
difference  in  the  overall  behavior  between  the  gas  and  the  particles. 

In  Figures  31,  36  and  41,  v'  and  v’  are  the  particle  relative  velocity 

F1  p2 

components  in  the  x-e  surface  at  the  inlet  and  the  exit  from  the  rotor, 
respectively,  while  and  repz*esent  their  corresponding  absolute 

values.  The  particle  radial  velocity  components  are  written  next  to 
the  corresponding  velocity  diagrams. 

The  effect  of  particle  material  density,  pp,  on  the  dynamic 
behavior  of  the  solid  particles  is  shown  in  Figures  42.  through  53. 

The  dynamic  behavior  of  the  three  typical  particles  is  shown  in  three 
groups  of  four  figures  each  for  pp  =  34,  68.7  and  151  ib/ft3  for  a 
constant  dp  =  40  microns  and  C  /C  =0.3.  These  particle  material 

J?i  9^ 

densities  approximately  correspond  to  the  densities  of  coal,  sand,  and 
silicon,  which  are  the  materials  that  a  practical  gas  turbine  flow  may 
contain  as  solid  suspension.  It  may  be  shown  from  the  above  figures 
that  a  decrease  in  Pp  will  lessen  the  deviation  between  the  particle 
trajectory  and  the  gas  streamlines  (Figures  42,  46  and  50),  decrease 
their  radial  displacements  (Figures  44,  48  and  52),  and  increase  their 
initial  acceleration  (Figures  45,  49  and  53).  Since  particxes  with 
smaller  pp  tend  to  escape  collision  with  the  rotor  blades,  they  ex¬ 
perience  fewer  abrupt  changes  in  their  velocity  as  they  pass  through 
the  rotor  (Figure  53) . 

3 

For  particles  with  constant  dp  =  40  microns  and  pp  =  68.7  Ib/ft  , 
the  initial  nondimens ion al  absolute  particle  velocity  C  /C  is  varied, 

pi  ^i 

taking  the  values  0.15,  0.3  and  0.6  to  determine  the  effect  of  C  /C  on 

Pi 
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the  particle  dynamic  behavior.  The  dynamic  behaviors  of  the  three 
typical  particles  are  plotted  in  Figures  54  through  65.  These  figures 
show  that  the  particle  initial  speed  have  a  small  effect  on  their 
trajectories,  especially  in  the  stator.  The  particles  are  found  to 
reach  the  same  absolute  velocity  after  traveling  relatively  short 
distances  irrespective  of  their  initially  different  value  of  C  /Cq 

pi  yi 

(Figure  57) .  The  effect  of  the  initial  gas  velocity  on  the  dynamic 
behavior  of  solid  particles  may  be  determined  from  the  figures  for 
different  C  /C  ,  bearing  in  mind  that  for  the  same  inlet  particle 
pi  "i 

velocity,  higher  gas  inlet  velocity  would  mean  lower  C  /C  ,  and  vice 
versa.  pi  ^i 


The  particle  mean  diameter  is  the  parameter_that  has  the  greatest 
effect  on  the  particles  dynamic  behavior,  while  pp  and  C  /C  have 

r  Pi  9i 


smaller  effects.  The  same  conclusion  could  be  arrived  at  by  looking 
at  the  equations  of  motion  of  the  particles,  since  d  appears  raised 
to  the  second  power.  p 


19 


CONCLUSION 


The  dynamic  behavior  of  solid  particles  entrained  by  the  compressible 
gas  flow  in  a  stationary  or  rotating  cascade  of  a  turbine  is  determined. 
Impacts  of  the  solid  particles  with  the  blades  and  the  casing  are  con¬ 
sidered.  The  equations  of  motion  of  the  solid  particles  are  solved  in 
the  three  dimensional  space.  The  drag  forces  on  the  particles  are 
calculated  using  drag  formulas  that  fit  the  drag  curve  of  a  spherical 
particle  over  a  wide  range  of  the  Reynolds  numbers.  The  compressible 
gas  flow  properties  are  computed  by  solving  numerically  the  gas  equations 
of  motion  in  a  blade  to  blade  surface  of  revolution  of  a  rotating  cas¬ 
cade.  Experimental  investigation  is  made  to  study  the  impact  and  re¬ 
bound  phenomenon  of  the  particles  from  the  walls,  to  determine  formulas 
for  the  restitution  ratio,  and  rebound  to  incidence  angle  ratio.  These 
formulas  define  the  particle  conditions  after  collision  which  are  then 
used  to  continue  the  solution  of  the  equations  of  motion  for  the  particles. 
These  formulas  have  to  be  determined  experimentally  for  every  particle- 
target  material  combination  once  the  turbomachine  operating  condition 
is  known. 

The  study  showed  that  in  general,  solid  particle  paths  are  deviated 
from  gas  streamlines.  This  deviation  increases  with  increased  particle 
mean  diameter,  material  density,  particle  initial  velocity  or  decreased 
gas  initial  velocity.  The  particle  mean  diameter  has  the  greatest 
effect  on  the  dynamic  behavior  of  the  particles,  while  the  particle 
material  density  has  a  lesser  effect  and  the  particle  initial  velocity 
has  the  least  effect. 

The  turbine  stator  and  rotor  blade  leading  edges  and  pressure  sides 
will  be  eroded  by  particles  of  all  sizes.  The  rear  part  of  the  stator 
suction  side  will  be  eroded  by  larger  particles  after  their  re*-\.rn  to 
the  stator  due  to  their  collision  with  the  rotor  olades.  The  turbine 
rotor  is  expected  to  suffer  more  severe  erosion  than  the  stator  due  to 
the  higher  velocities  of  the  particles  in  the  rotor. 

Turbines  may  be  designed  to  both  minimize  erosion  as  well  as 
optimizing  aerodynamic  characteristics  using  the  results  of  this  in¬ 
vestigation.  Means  may  be  introduced  to  collect  or  deviate  some  of 
the  particles,  expecially  those  that  contribute  most  to  erosion,  from 
the  blades  and  hence,  reduce  erosion  damage. 
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FIGURE  12  TURBINE  STAGE  DIMENSIONS  FIGURE  13  COMBINED  GAS  VELOCITY  DIAGRAM. 
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FIGURE  14  MESH  POINTS  FOR  TURBINE  BLADING 
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FIGURE  17  AXIAL  AND  RADIAL  COMPONENTS 
OF  PARTICLE  TRAJECTORIES 


FIGURE  18  PARTICLE  NONDIMENSIONAL  ABSOLUTE 
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FIGURE  30  PARTICLE  NONDIMEN  SION  AL  ABSOLUTE  FIGURE  31  PARTICLE  VELOCITY  DIAGRAM 
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FIGURE  33  AXIAL  AND  TANGENTIAL  COMPONENTS 
OF  PARTICLE  TRAJECTORIES  RELATIVE 
TO  THE  ROTOR  BLADES 
(  EFFECT  OF  dp ) 
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FIGURE  35  PARTICLE  NONDIMENSiONAL  ABSOLUTE  -  FIGURE  36  PARTICLE  VELOCITY  DIAGRAM 
VELOCITIES  (EFFECT  OF  dD)  • 
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FIGURE  40  PARTICLE  NONDIMENSIONAL  ABSOLUTE  FIGURE  41  PARTICLE  VELOCITY  DIAGRAM 
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FIGURE  45  PARTICLE  NONDIMENSIONAL  ABSOUUTE 
VELOCITIES  (EFFECT  OF  Pp) 
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FIGURE  60  AXIAL  AND  RADIAL  COMPONENTS 
OF  PARTICLE  TRAJECTORIES 


FIGURE  61  PARTICLE  N  ONDIMENSiONAL  ABSOLUTE 
VELOCITIES  (EFFECT  OF  Cp./Cg.) 
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